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Abstract 

We study dynamics of a ball moving in gravitational field and colliding 
with a moving table. The motion of the limiter is assumed as periodic 
with piecewise constant velocity - it is assumed that the table moves 
up with a constant velocity and then moves down with another constant 
velocity. The Poincare map, describing evolution from an impact to the 
next impact, is derived and scenarios of transition to chaotic dynamics 
are investigated analytically. 

1 Introduction 

Vibro-impacting systems belong to a very interesting and important class of 
nonsmooth and nonlinear dynamical systems [1] [2J 131 H] with important tech- 
nological applications [SI [51 [Zl IH] • Dynamics of such systems can be extremely 
complicated due to velocity discontinuity arising upon impacts. A very charac- 
teristic feature of such systems is the presence of nonstandard bifurcations such 
as border-collisions and grazing impacts which often lead to complex chaotic 
motions. 

The Poincare map, describing evolution from an impact to the next impact, 
is a natural tool to study vibro-impacting systems. The main difhculty with 
investigating impacting systems is in finding instant of the next impact what 
typically involves solving a nonlinear equation. However, the problem can be 
simplified in the case of a bouncing ball dynamics assuming a special motion 
of the limiter. In the present paper we investigate motion of a material point 
in a gravitational field colliding with a limiter moving with piecewise constant 
velocity. This class of models has been extensively studied, see [9] and references 
therein. As a motivation that inspired this work, we mention study of physics 
and transport of granular matter [6]. A similar model has been also used to 
describe the motion of railway bogies [7]. Therefore it can be expected that 
some of the present results may cast light on the dynamics in such systems. 
On the other hand, simple motion of the limiter makes analytical explorations 
possible, cf. our preliminary report |10| . 
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The paper is organized as follows. In Section 2 a one dimensional dynamics of 
a ball moving in a gravitational field and colliding with a table is considered and 
Poincare map is described for piecewise linear motion of the table. In Section 3 
transition to chaotic dynamics from periodic motion is described. The nature of 
mixing leading to chaotic dynamics is described in Section 4. In Sections 5 and 
6 homoclinic structures responsible for mixing are determined and computed. 
Finally, we discuss our results in the last Section. 



2 Bouncing ball: a simple motion of the table 

We consider a motion of a small ball moving vertically in a gravitational field 
and colliding with a moving table, representing unilateral constraints. The ball 
is treated as a material point while the limiter's mass is assumed so large that 
its motion is not affected at impacts. A motion of the ball between impacts is 
described by the Newton's law of motion: 

mx = —mg, (1) 

where x = dx/dt and motion of the limiter is: 

y^y{t), (2) 

with a known function y. We shall also assume that y is a continuous function 
of time. Impacts are modeled as follows: 

x{Ti) = 2/(Ti), (3) 

x{Tt)-y{n) - ^R{x{Tr)-y{n)), (4) 

where duration of an impact is neglected with respect to time of motion between 
impacts. In Eqs. (jSj, (|4]) stands for time of the i-th impact while i~, i^are 
left-sided and right-sided limits of Xi (t) for t ^ Ti, respectively, and R is the 
coefficient of restitution, < i? < 1 5 . 

Solving Eq. ([T]) and applying impact conditions (jS]), ([4]) we derive the 
Poincare map 11 : 

V,+i = -RV, + 2RA,+i +j{l + R)Y (T,+i) , ^ ' 

where Aj+i = T^+i — Tj. The limiter's motion has been typically assumed 
in form Y{T) — sin(r), cf. [T^] and references therein. This choice leads to 
serious difficulties in solving the first of Eqs. ([5]) for Ti+i, thus making analytical 
investigations of dynamics hardly possible. Accordingly, we have decided to 
simplify the limiter's periodic motion to make ^ solvable. Let us thus assume 
that the table moves up with a finite constant velocity 7Y1 and then goes down 
with a finite constant velocity 7I2 !T2]- Therefore, displacement of the table is 
the following periodic function of time: 

^ T^kiT-lTD + jh^ {h<T-lT\) 
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with Yi — l/hjYz = —1/ {1 — h), < h < 1, where [a;J is the floor function - the 
largest integer less than or equal to x. Our model consists thus of equations (O, 
^ with control parameters R, 7, h. Since the period of motion of the limiter 
is equal to one, the map ([S]) is invariant under the translation Ti Ti + 1. 
Accordingly, all impact times Ti can be reduced to the unit interval [0, 1]. 



3 From periodic dynamics to chaotic motion 



In our recent article periodic solutions of Eqs. (O, dH) have been investigated 
[T^ . Dynamics of Eqs. (0, (O becomes complicated when some impacts occur 
in time interval (0, h), and some in {h, 1). 

In the case of two impacts per period a 2 - cycle, Ti e (0, h), T2 G {h, 1), 
is stable. Next, for increasing values of 7, period doubling takes place and 2^ - 
cycle with impacts < T*i,T*3 < h, h < r*2,T'*4 < 1 is formed. Then, upon 
further increase of 7, the period doubling scenario ends abruptly when T^2 = 1- 
For R = 0.85, h = 0.2623 this happens for 7^^ = 0.09058194712. This critical 
transition is investigated in Sections 4, 5. Equations for dynamics after the first 
period doubling have the following form: 



7( 



= IT, - (T2 



i-h - h^^ v-^^ ^1) 
V2 = -RVi + 2R{T2-Ti)- 

7(T3-1) _ 7(-T2 + l) (rr, 

V3 = -RV2 + 2R (T3 - T2) + 
7(-r4+i) - 2I3. - (T - T ^2 

1-h ~ h -^3; 

Vi = -RV3 + 2R{Ti-n) + 



+ {T2- 

7(l + fl) 



Ti)Vi 



1-h 

T2f + (Ta 



T2)V2 



7(l + fl) 
h 

+ {T4- 
7(1 + -R) 



^3)^3 



(7) 



7(l + fl) 
h 



7^5 _ 7^5 (rp rp \2 

V5 = -RV4 + 2R{T5-T4) 

^ V5 = Vi 

In Fig. 1 transition to chaos is shown. The initial dynamical state with 
two impacts per period bifurcates at 7 = 7^^ = 0.090099 (this value can be 

computed analytically, see [l2]). For 7 7^"* time of the second impact tends 

(2) 

to 1 and this mode of dynamics is impossible for 7 > jcr ■ It turns out that for 
(2) 

7 > 7cr there are two attractors: a noisy, probably chaotic, attractor coexisting 
with a 7-cycle which appears just before the transition. At 7 = 0.09060 the 
noisy attractor disappears and is substituted by a more irregular attractor, see 
Figs. 2,3. 

Full circles indicate positions of small clouds of points. We have studied these 
potentially chaotic attractors in detail. First of all, we have checked numerically 
that the attractors are non-periodic. Indeed, computations show that after 10^ 
iterations the points generated by the map (O, ([6]) stay on the corresponding 
attractor and do not repeat. We have also computed Lyapunov exponents for 
both attractors. In the case of the attractor shown in Fig. 2 the Lyapunov 
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Figure 1: Bifurcation diagram: R = 0.85, h = 0.2623, 7 G [0.09, 0.090605]. 



exponent is Ai = 0.5 while for the attractor in Fig. 3 A2 = 0.8. It foUows that 

in both cases dynamics is chaotic and is more mixing in the second case. The 
mechanism of mixing is explained in the next Section. 



"V 










■ 1 





Figure 2: Chaotic attractor, R = 0.85, h = 0.2623, 7 = 0.09059. 
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Figure 3: Chaotic attractor, R = 0.85, h = 0.2623, 7 = 0.09060. 

4 Mechanism of mixing 

Mixing can arise due to corner events [T] when impacts occur at points where 
motion of the hmiter loses smoothness at time instances Tcl^ — Tcr^ — 1. 
Let us investigate the second possibihty more closely. In Fig. 4 the stable 2^ - 
cycle with four impacts per two periods: < T*i, T*3 < h, h < T^2, T*4 < 1 and 
unstable 2 - cycle are shown schematically. 




Figure 4: Stable 2^ - cycle (circles) and the unstable 2 - cycle (boxes). 

(2) 

For increasing value of the control parameter 7 we have T,2 ^ Tcr = 1, 
see Fig. 1. The map ([5|), ([6]) is invariant under translation Ti — > T^+i and 
thus the phase space is topologically equivalent to the cylinder and hence we 

have to glue the end points of the time interval [0, 1] obtaining thus a circle. 

(2) (2) 
Therefore, a small neig hborhood of TcV^ = 1 is a union of two sets, Ocr = 

{T : 1 - ei < T < l}U{r : < T < 62}, where ei, 62 are smah and positive. 
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see Fig. 5. 
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Figure 5: Unstable 2^ - cycle (circles) and the unstable 2 - cycle (boxes). 

Now, let Ti e 0il\ It follows from Eqs. ©, © that time of the next 
impact, Ti+i, as well as the corresponding post-impact velocity, Vi^i, depend 
discontinuously on Ti. In other words, we get different solutions, Ti^i, 14+ 1, 
depending on whether 1 — ei < < 1 or < < £2- It follows that mixing 
will necessarily be present if a trajectory recurrently visits the interval Ocr ■ We 
shall study this possibility in the next Section. 

5 Origin of the homoclinic structure 

Let us consider critical case: 7 = Jc^\ r*2 = Tc? — 1- This is described 
by Eq. ([71) with four impacts per two periods such that < T*i,T*3 < ft,, 



where we have substituted T2 = 1 and the periodicity conditions: Ts Ti + 1, 
V5 = V\. Solutions of Eq. ([8]) depend on roots of the following algebraic 
equation: 



' = ^ri-(i-ri)V(i-Ti)T/i 

F2 = -RVr + 2i? (1 - Ti) - liii^ 
^ = - (T3 - 1) + ^2 



< 




(8) 



AgX^ + A^X^ + A^X^ + ^3^3 + ^2^2 + A^X + = 0, 



(9) 
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where X = Ti with coefScients Ai depending on R, h in a complicated way. 
The coefficients are Usted below: 

Ae = {R+lf{R^ + lf, 



^5 = as + hh: 

as = 2{R+1) 

65 = -4i? {R^ -R- 2) {R + if [R^ + 1) 



-48i?5 _ 36i?4 + 23i?-^ - 9i?2 - Hi? - 3 



— /l/'D_L1^^ — ii-fx — 'ii- 



A4 = 04 + 64/1 + 04/1^: 

a4 = 23 + 74i? - 6i?2 - 158i?,^ + 121i?/' + 348i?^ + 316i?'' 

+28i?^ + 73i?* + 250i?^ + 170i?i° + 34i?" + IR^'^, 

, _ / 17i?io + llii9-20ii8-49E7 + 5EM 

b4 - 4(K+lj(^ -19i?5-41i?4-55i?3 + 10i?2-3 j' 

C4 = A{R^ + l){R+lf{R^ + l)\ 
^3 = as + '>3/i + C3/i^: 

i?" - _ 4ii?9 _ 55^8 _^ 287?,^ - 32i?6 

-48i?'"^ - 32i?4 + 3i?3 + 43i?2 _ 7i? _ 9 
63 = 16 - 72i? - 160i?2 + 352i?3 + 5207?^ + 336i?5 + 32i?6 

+416i?^ + hl2R^ + 184i?^ - 128i?^° - 64i?" - 247?^^ 
C3 = -8i?(i?+l)(i?'' + l) (i?'5+9i?^ + 2i?*-i?2 + 15ii-2), 

A2 = a2 + 62/1 + C2ft.^ + (^2/1^: 

02 = 31 - 2i? - 118i?2 _ i78i?3 _^ 9i?* + 188i?5 + 252i?^ 
-ISOi?^ + 113i?^ + 358i?^ + 154i?^° + 6i?" + 7i^^^ 

62 = -8 + 192i? + 120i?2 - 168i?^ - AIQR^ - 176i?^ - 192i?^ 

-632i?^ - 208i?^ + 72i?^° + 40i?" - 32i?^^ 

6i?/^ + 10i?^ + lli?^-8i?5 
-'iR'^ + 14i?-* + 21i?2 - 8i? - 1 

d2 = 16i?(i?-l)(i^+l)(i?'' + l)^ 

Ai = ai + 61/1 + ci/i^ + dih^: 

3i?i° + 4i?9 - 3i?8 + 92i?^ - 132i?6 
+80i?5 - 16i?4 + 28i?3 - 47i?2 + 4i? + 3 

61 = -32 - 64i? - 180i?^ + 56i?^ + 204i?^ + 1127?^ - 280i?'' 

+128i?'' + 392i?^ + 80i?^ - 116i?^" + 72i?" + 127?^^ 
ci = 8 (i?" + 1) {R^ - 12i?^ + R^ + 8R^ + R^- 20R^ -R'^ + 6), 

di = -16{R-l){R^ + R'^ + R-l){R'^ + lY , 



C2 = 8{R^ + 1) 



ai = -2{R+lf 
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^0 = ao + bah + coh'^ + doh^: 



ao = 

bo = 

Co = 

do = 



{R - 1) {R^ + 4i?2 - 

-4(i?-l) 

-4 {R - 1) + 1) 
16R{R-lf (i^-^ + l 



4i? + 3) 



i?^ + 7R^ + 9R^ - ISi?-* 
+7R^ + 9i?2 _ _ 3 

9i?4 + 23i?3 + 24i?2 + 8i? - 3 

+3i?3 - 7i?2 - 13i? + 1 
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Let us stress here that acceptable solution for the time of the first impact 
must fulfill consistency condition Ti = X € (0, h). Now it follows that necessary 
conditions for existence of this solution can be formulated. Indeed, the condition 
Ao = f {R, h) ^ guarantees existence of solution Ti = X ~ 0. Furthermore, 
after change of variable X = X + h the equation ^ is written as AqX^ + . . . + 
AiX + Ao =0 and the condition Aq = AqH^ + A^h^ +^4/1'' + A^h^ + ^2/1^ + 
Aih + Ao — guarantees existence of the solution X = and hence existence 
of solution Ti = X = h. Region of acceptable values of parameters R, h is 
shown in Fig. 6 - it is placed between thin solid lines (which correspond to the 
condition Ti = h) and below medium solid line (the condition Ti =0). 




Figure 6: Acceptable parameter values (shaded region). 

The solution of Eq. ([8]) is unstable and leads to a homoclinic-type orbit 
and thus can be referred to as the homoclinic cycle (see [T51 for the definition 
of a homoclinic point and a homoclinic orbit). Indeed, for the initial condition 
1 - £1 < T2 < 1 the orbit is attracted by the 2^ - cycle r*i,r*2 = l,r*3,r*4, 
while for < T2 < €2 the fixed point T»2 = 1 is repelling, but the orbit returns 
eventually to the 2^ - cycle (provided that a coexisting attractor does not capture 
the trajectory). We shall compute the repelling branch in the next Section. 
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6 Computing the homoclinic orbit 



To analyse structure of the chaotic attractor shown in Fig. 2 we have solved Eq. 
([8]) for i? = 0.85, h = 0.2623, computing thus critical value of the parameter 7 
and the homoclinic cycle where (r,2, ^^2) is the homoclinic point: 

^ 0.090 581947119, 

T,i = 0.122 922181823, Ki ==0.828 678 917871, 

T,2 = l", 142 = 0.559 494 302 251, (10) 
T,3 = 0. 214 157 103 932 , = 0.527 370 736 661 , 

T,4 = 0.827 577 663 824 , V^i = 0.367 388 917 196 . 

This solution is attracting for initial condition T2 £ (1 — ei,l) and V2 = 
y*2- The sequence {Ti,Vi} starting from such initial condition belongs to the 
attracting branch of the homoclinic orbit. The repelling branch is obtained in 
the following way. We start from Ti = T^i, V\ =V*\. T2 and V2 are computed 
from the following equations: 

^(T2-i) = Jr*i-(T2-T,i)V(r2-T,i)Ki (11) 

V2 = -RV,i + 2R{T2-T,i) + [l + R)l (12) 

h 



The solution of the first equation is of course T2 = 1. We assume now in (|12l) 
that at the impact the table is just about going up with velocity Y {T2) = 1/h 
rather than it has just finished going down with velocity Y (T2) = —1/ {1 — h) 
(therefore this equation differs from the second of equations in ([8|)). We thus 
compute from Eq. ([12]) , for i? = 0.85, h = 0.2623 and 7 = 7^^^ = 0.090 581 947 119, 



that V2 = 1. 425 529 027 600 (let us stress again that using the second of Eqs. ^ 
we get = 0.559 494 302 251) . Due to symmetry of the dynamics Ti^Ti + 1 
the first point of the repelling branch of the homoclinic orbit can be assumed 
as Tf = 0, = 1. 425 529 027 600 . 

We have thus computed numerically the repelling branch of the homoclinic 
orbit starting from the initial condition T-^, V;^ {R = 0.85, h = 0.2623, 7 = 
-fil'' = 0.090 581947119). 

We have shown in Fig. 6 first twenty six points (full circles) of the repelling 
branch of the homoclinic orbit starting from the point (T/^, V^) - the outermost 
full circle in the figure, lying on the vertical axis. These points agree very well 
with positions of twenty six clouds of points belonging to the chaotic attractor 
shown in Fig. 2, computed for 7 = 0.09059. The next points (dots) of the 
repelling branch of the homoclinic trajectory enter four connected parts placed 
as in Fig. 2 and tend, as an attracting branch, to the homoclinic cycle (larger 
open circles) containing the homoclinic point (T*2,142) - the outermost open 
circle. 

This homoclinic structure is preserved in the interval 7cr'' = 0.0905819471 . . . < 
7 < 0.09060, R — 0.85, h = 0.2623 until this attractor is substituted by a new 
one due to crisis (the unstable cycle collides with one of clouds of points belong- 
ing to the attractor), see the bifurcation diagram. Fig. 1, and Figs. 2,3. 
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Figure 7: The homoclinic orbit. 

7 Summary and discussion 

We have found a generic scenario of transition to chaos for dynamics of a baU 
moving verticahy in gravitational field and colliding with a table moving verti- 
cally with piecewise constant velocity. 

According to this scenario a periodic and stable solution is destroyed via a 
corner bifurcation [1] in a corner event, T^i = Tcr^ = h or T^i — Tc? = 1. In 
the present paper the solution, defined analytically by T*2 = 1 in Eq. (O, is a 
homoclinic-type orbit and leads to mixing and hence to chaotic dynamics. This 
homoclinic-type orbit is untypical in the sense that it is not a saddle structure 

(2) 

but its origin is related to discontinuous dynamics in the neighborhood of TcV . 
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